Catalytic antibody response to the Staphylococcus aureus virulence factor Efb

Staphylococcus aureus is a globally prevalent pathogen that can cause a wide variety of acute and chronic diseases in both adults and children, in both immune susceptible populations and healthy individuals. Its ability to cause persistent infections has been linked to multiple immune evasion strategies, including Efb-mediated complement inhibition. As new multi-drug-resistant strains emerge, therapeutic alternatives to traditional antibiotics must be developed. These experiments assessed the ability of healthy patient immunoglobulin to cleave Efb and disable the complement-inhibitory properties of Efb in vitro. Levels of immunoglobulin-mediated Efb catalysis varied both between immunoglobulin isoform/isotype and between individuals. Serum IgG showed the strongest catalytic activity of the immunoglobulin isotypes tested. Additionally, IgG hydrolyzed the virulence factor in a way that enabled only minimal binding to the complement component C3b, effectively blocking Efb-mediated inhibition of complement lysis. Salivary IgA and serum IgM did not block Efb-mediated inhibition of complement. Catalytic IgG selectively cleaved Efb and showed no cleavage of a variety of other proteins tested. Catalytic activity of IgG was inhibited by serine protease inhibitors, but not by other protease inhibitors, suggesting a serine-protease mechanism of catalysis. It is proposed that varying concentrations and activity levels of catalytic IgG between healthy individuals and those with current or recurrent S. aureus infections in both adult and pediatric populations be studied in order to assess the potential effectiveness of passive immunization therapy with catalytic monoclonal IgG. ^

Antibody response to Panton-Valentine Leukocidin, a Staphylococcus aureus virulence factor, in pediatric patients from Texas Children's Hospital in Houston, Texas

Staphylococcus aureus is an important human pathogen of global health significance, whose frequency is increasing and whose persistence and versatility allow it to remain established in communities worldwide. An observed significant increase in infections, particularly in children with no predisposing risk factors or medical conditions, led to an investigation into pediatric humoral immune response to Panton-Valentine Leukocidin (PVL) and to other antigens expressed by S. aureus that represent the important classes of virulence activities. Patients who were diagnosed with staphylococcal infections were enrolled (n=60), and serum samples collected at the time of admission were analyzed using ELISA and Western blot to screen for immune response to the panel of recombinant proteins. The dominant circulating immunoglobulin titers in this pediatric population were primarily IgG, were specific, and were directed against LukF and LukS, while suppression of other important virulence factors in the presence of PVL was suggested. Patients with invasive infections (osteomyelitis, pneumonia or myositis) had higher titers against LukF and LukS compared to patients with non-invasive infections (abscesses, cellulitis or lymphadenitis). In patients with osteomyelitis, antibody responses to LukF and LukS were higher than antibody responses to any other virulence factor examined. This description of immune response to selected virulence factors of S. aureus caused by isolates of the USA300 lineage in children is novel. Antibody titers also correlated with markers of inflammation. The significance of these correlations remains to be understood.^

Control of the master virulence regulatory gene atxA in Bacillus anthracis

Transcription of the Bacillus anthracis structural genes for the anthrax toxin proteins and biosynthetic operon for capsule are positively regulated by AtxA, a transcription regulator with unique properties. Consistent with the role of atxA in virulence factor expression, a B. anthracis atxA-null mutant is avirulent in a murine model for anthrax. In batch culture, multiple signals impact atxA transcript levels, and the timing and steady state level of atxA expression is critical for optimal toxin and capsule synthesis. Despite the apparent complex control of atxA transcription, only one trans-acting protein, the transition state regulator AbrB, has been demonstrated to directly interact with the atxA promoter. The AbrB-binding site has been described, but additional cis-acting control sequences have not been defined. Using transcriptional lacZ fusions, electrophoretic mobility shift assays, and Western blot analysis, the cis-acting elements and trans-acting factors involved in regulation of atxA in B. anthracis strains containing either both virulence plasmids, pXO1 and pXO2, or only one plasmid, pXO1, were studied. This work demonstrates that atxA transcription from the major start site P1 is dependent upon a consensus sequence for the housekeeping sigma factor SigA, and an A+T-rich upstream element (UP-element) for RNA polymerase (RNAP). In addition, the data show that a trans-acting protein(s) other than AbrB negatively impacts atxA transcription when it binds specifically to a 9-bp palindrome within atxA promoter sequences located downstream of P1. Mutation of the palindrome prevents binding of the trans-acting protein(s) and results in a corresponding increase in AtxA and anthrax toxin production in a strain- and culture-dependent manner. The identity of the trans-acting repressor protein(s) remains elusive; however, phenotypes associated with mutation of the repressor binding site have revealed that the trans-acting repressor protein(s) indirectly controls B. anthracis development. Mutation of the repressor binding site results in misregulation and overexpression of AtxA in conditions conducive for development, leading to a marked sporulation defect that is both atxA- and pXO2-61-dependent. pXO2-61 is homologous to the sensor domain of sporulation sensor histidine kinases and is proposed to titrate an activating signal away from the sporulation phosphorelay when overexpressed by AtxA. These results indicate that AtxA is not only a master virulence regulator, but also a modulator of proper B. anthracis development. Also demonstrated in this work is the impact of the developmental regulators AbrB, Spo0A, and SigH on atxA expression and anthrax toxin production in a genetically incomplete (pXO1+, pXO2-) and genetically complete (pXO1+, pXO2+) strain background. AtxA and anthrax toxin production resulting from deletion of the developmental regulators are strain-dependent suggesting that factors on pXO2 are involved in control of atxA. The only developmental deletion mutant that resulted in a prominent and consistent strain-independent increase in AtxA protein levels was an abrB-null mutant. As a result of increased AtxA levels, there is early and increased production of anthrax toxins in an abrB-null mutant. In addition, the abrB-null mutant exhibited an increase in virulence in a murine model for anthrax. In contrast, virulence of the atxA promoter mutant was unaffected in a murine model for anthrax despite the production of 5-fold more AtxA than the abrB-null mutant. These results imply that AtxA is not the only factor impacting pathogenesis in an abrB-null mutant. Overall, this work highlights the complex regulatory network that governs expression of atxA and provides an additional role for AtxA in B. anthracis development.

Evidence for VirB4-mediated dislocation of membrane-integrated VirB2 pilin during biogenesis of the Agrobacterium VirB/VirD4 type IV secretion system.

Agrobacterium VirB2 pilin is required for assembly of the VirB/VirD4 type IV secretion system (T4SS). The propilin is processed by signal sequence cleavage and covalent linkage of the N and C termini, and the cyclized pilin integrates into the inner membrane (IM) as a pool for assembly of the secretion channel and T pilus. Here, by use of the substituted cysteine accessibility method (SCAM), we defined the VirB2 IM topology and then identified distinct contributions of the T4SS ATPase subunits to the pilin structural organization. Labeling patterns of Cys-substituted pilins exposed to the membrane-impermeative, thiol-reactive reagent 3-(N-maleimidopropionyl)biocytin (MPB) supported a topology model in which two hydrophobic stretches comprise transmembrane domains, an intervening hydrophilic loop (residues 90 to 94) is cytoplasmic, and the hydrophilic N and C termini joined at residues 48 and 121 form a periplasmic loop. Interestingly, the VirB4 ATPase, but not a Walker A nucleoside triphosphate (NTP) binding motif mutant, induced (i) MPB labeling of Cys94, a residue that in the absence of the ATPase is located in the cytoplasmic loop, and (ii) release of pilin from the IM upon osmotic shock. These findings, coupled with evidence for VirB2-VirB4 complex formation by coimmunoprecipitation, support a model in which VirB4 functions as a dislocation motor to extract pilins from the IM during T4SS biogenesis. The VirB11 ATPase functioned together with VirB4 to induce a structural change in the pilin that was detectable by MPB labeling, suggestive of a role for VirB11 as a modulator of VirB4 dislocase activity.

Studies of subcellular localization and complex formation of the Agrobacterium tumefaciens transport ATPase VirB11

The VirB11 ATPase is an essential component of an Agrobacterium tumefaciens type IV bacterial secretion system that transfers oncogenic nucleoprotein complexes to susceptible plant cells. This dissertation investigates the subcellular localization and homo-oligomeric state of the VirB11 ATPase in order to provide insights about the assembly of the protein as a subunit of this membrane-associated transfer system. Subcellular fractionation studies and quantitative immunoblot analysis demonstrated that $\sim$30% of VirB11 partitioned as soluble protein and $\sim$70% was tightly associated with the bacterial cytoplasmic membrane. No differences were detected in VirB11 subcellular localization and membrane association in the presence or absence of other transport system components. Mutations in virB11 affecting protein function were mapped near the amino terminus, just upstream of a region encoding a Walker 'A' nucleotide-binding site, and within the Walker 'A' motif partitioned almost exclusively with the cytoplasmic membrane, suggesting that an activity associated with nucleotide binding could modulate the affinity of VirB11 for the cytoplasmic membrane. Merodiploid analysis of VirB11 mutant and truncation derivatives provided strong evidence that VirB11 functions as a homo- or heteromultimer and that the C-terminal half of VirB11 contains a protein interaction domain. A combination of biochemical and molecular genetic approaches suggested that VirB11 and the green fluorescence protein (GFP) formed a mixed multimer as demonstrated by immunoprecipitation experiments with anti-GFP antibodies. Second, a hybrid protein composed of VirB11 fused to the N-terminal DNA-binding domain of bacteriophage $\lambda$ cI repressor conferred immunity to $\lambda$ superinfection, demonstrating that VirB11 self-association promotes dimerization of the chimeric repressor. A conserved Walker 'A' motif, though required for VirB11 function in T-complex export, was not necessary for VirB11 self-association. Sequences in both the N- and the C-terminal halves of the protein were found to contribute to self-association of the full length protein. Chemical cross-linking experiments with His$\sb6$ tagged VirB11 suggested that VirB11 probably assembles into a higher order homo-oligomeric complex. ^

AGROBACTERIUM VIRB10 CONTRIBUTIONS TO TYPE IV SUBSTRATE SECRETION, T-PILUS BIOGENESIS, AND OUTER MEMBRANE PORE FORMATION

The VirB/D4 type IV secretion system (T4SS) of Agrobacterium tumefaciens functions to transfer substrates to infected plant cells through assembly of a translocation channel and a surface structure termed a T-pilus. This thesis is focused on identifying contributions of VirB10 to substrate transfer and T-pilus formation through a mutational analysis. VirB10 is a bitopic protein with several domains, including a: (i) cytoplasmic N-terminus, (ii) single transmembrane (TM) α-helix, (iii) proline-rich region (PRR), and (iv) large C-terminal modified β-barrel. I introduced cysteine insertion and substitution mutations throughout the length of VirB10 in order to: (i) test a predicted transmembrane topology, (ii) identify residues/domains contributing to VirB10 stability, oligomerization, and function, and (iii) monitor structural changes accompanying energy activation or substrate translocation. These studies were aided by recent structural resolution of a periplasmic domain of a VirB10 homolog and a ‘core’ complex composed of homologs of VirB10 and two outer membrane associated subunits, VirB7 and VirB9. By use of the substituted cysteine accessibility method (SCAM), I confirmed the bitopic topology of VirB10. Through phenotypic studies of Ala-Cys insertion mutations, I identified “uncoupling” mutations in the TM and β-barrel domains that blocked T-pilus assembly but permitted substrate transfer. I showed that cysteine replacements in the C-terminal periplasmic domain yielded a variety of phenotypes in relation to protein accumulation, oligomerization, substrate transfer, and T-pilus formation. By SCAM, I also gained further evidence that VirB10 adopts different structural states during machine biogenesis. Finally, I showed that VirB10 supports substrate transfer even when its TM domain is extensively mutagenized or substituted with heterologous TM domains. By contrast, specific residues most probably involved in oligomerization of the TM domain are required for biogenesis of the T-pilus.

Characterization of the Agrobacterium tumefaciens VirB2 pilin of the VirB/D4 Type IV Secretion System

The Agrobacterium tumefaciens VirB/D4 type IV secretion system (T4SS) delivers oncogenic T-DNA and effector proteins to susceptible plant cells. This leads to the formation of tumors termed Crown Galls. The VirB/D4 T4SS is comprised of 12 subunits (VirB1 to VirB11 and VirD4), which assemble to form two structures, a secretion channel spanning the cell envelope and a T-pilus extending from the cell surface. In A. tumefaciens, the VirB2 pilin subunit is required for assembly of the secretion channel and is the main subunit of the T-pilus. The focus of this thesis is to define key reactions associated with the T4SS biogenesis pathway involving the VirB2 pilin. Topology studies demonstrated that VirB2 integrates into the inner membrane with two transmembrane regions, a small cytoplasmic loop, and a long periplasmic loop comprised of covalently linked N and C termini. VirB2 was shown by the substituted cysteine accessibility method (SCAM) to adopt distinct structural states when integrated into the inner membrane and when assembled as a component of the secretion channel and the T-pilus. The VirB4 and VirB11 ATPases were shown by SCAM to modulate the structural state of membrane-integrated VirB2 pilin, and evidence was also obtained that VirB4 mediates extraction of pilin from the membrane. A model that VirB4 functions as a pilin dislocase by an energy-dependent mechanism was further supported by coimmunoprecipitation and osmotic shock studies. Mutational studies identified two regions of VirB10, an N-terminal transmembrane domain and an outer membrane-associated domain termed the antennae projection, that contribute selectively to T-pilus biogenesis. Lastly, characterization of a VirB10 mutant that confers a ‘leaky’ channel phenotype further highlighted the role of VirB10 in gating substrate translocation across the outer membrane as well as T-pilus biogenesis. Results of my studies support a working model in which the VirB4 ATPase catalyzes dislocation of membrane-integrated pilin, and distinct domains of VirB10 coordinate pilin incorporation into the secretion channel and the extracellular T-pilus.

NMR structure of a complex between the VirB9/VirB7 interaction domains of the pKM101 type IV secretion system.

Type IV secretion (T4S) systems translocate DNA and protein effectors through the double membrane of Gram-negative bacteria. The paradigmatic T4S system in Agrobacterium tumefaciens is assembled from 11 VirB subunits and VirD4. Two subunits, VirB9 and VirB7, form an important stabilizing complex in the outer membrane. We describe here the NMR structure of a complex between the C-terminal domain of the VirB9 homolog TraO (TraO(CT)), bound to VirB7-like TraN from plasmid pKM101. TraO(CT) forms a beta-sandwich around which TraN winds. Structure-based mutations in VirB7 and VirB9 of A. tumefaciens show that the heterodimer interface is conserved. Opposite this interface, the TraO structure shows a protruding three-stranded beta-appendage, and here, we supply evidence that the corresponding region of VirB9 of A. tumefaciens inserts in the membrane and protrudes extracellularly. This complex structure elucidates the molecular basis for the interaction between two essential components of a T4S system.

Diffusely adherent Escherichia coli (DAEC) as a cause of acute diarrhea

Acute diarrhea is the most common medical problem in the developing countries. Infectious agents are responsible for a majority of cases of acute diarrhea. Knowing the cause of acute diarrhea is important to developing plans for disease prevention, control and therapy. Acute diarrhea is caused by many viruses, bacteria and parasites. ^ Travelers to developing countries of the world commonly develop diarrhea as a result of eating contaminated food or drinking contaminated water. About 30-50% of travelers who travel from industrialized countries like United States to the developing countries are at risk of developing diarrhea. High risk areas for travelers' diarrhea are Mexico, Latin America and Southeast Asia. Public restaurants are the common sites for exposure to this type of food-borne infectious disease in travelers. Food becomes contaminated when they are handled by people with fecal content on their hands. ^ The importance of Diffusely Adherent Escherichia Coli (DAEC) in travelers to these areas has not been well studied. Some of the studies looking at DAEC have shown the organism to be present in children without symptoms. Other studies have shown a relationship between DAEC infection and presence of symptoms. I have selected this topic because the patho-physiological processes in DAEC infection that allow intestinal and extra-intestinal infections to develop are not fully understood. DAEC related acute diarrhea is a relatively new topic of public health significance. There is a limited number of studies regarding the virulence and pathogenic mechanisms of DAEC. The presumed virulence factor of the organism is diffuse attachment to the intestinal lining of the infected host. However more research needs to be done to identify the pathogenic mechanisms and virulence factors associated with DAEC infection for better treatment planning and diarrhea prevention. ^

Regulations of SOX9 activity during chondrogenesis

Sox9 is a transcription factor required for chondrocyte differentiation and cartilage formation. In an effort to identify SOX9 interacting protein(s), we screened a chondrocyte cDNA library with a modified yeast two-hybrid method, Son of Sevenless (SOS) recruitment system (SRS). The catalytic subunit of cyclic AMP-dependent protein kinase A (PKA-Cα) and a new long form of c-Maf transcription factor (Lc-Maf) were found to interact specifically with SOX9. We showed here that two PKA phosphorylation consensus sites of SOX9 could be phosphorylated by PKA in vitro as well as in vivo. PKA phosphorylation of SOX9 increases its DNA binding and transcriptional activities on a Col2a1 chondrocyte-specific enhancer. Mutations of these two PKA phosphorylation sites markedly decreased the activation of SOX9 by PKA. ^ To test whether parathyroid hormone-related peptide (PTHrP) signaling results in SOX9 phosphorylation, we generated a phosphospecific antibody that specifically recognizes SOX9 that is phosphorylated at serine 181 (S 181) one of the two consensus PKA phosphorylation sites. Addition of PTHrP to COS7 cells cotransfected with SOX9 and PTH/PTHrP receptor strongly increased phosphorylation of SOX9 at S181; this phosphorylation was blocked by a PKA-specific inhibitor. In similar experiments we showed that PTHrP increased the activity of a SOX9-dependent Col2a1 enhancer. This increase in activity was abolished when a SOX9 mutant was used containing serine-to-alanine substitution in the two consensus PKA phosphorylation sites of SOX9. Using our phosphospecific SOX9 antibody we showed by immunohistochemistry of mouse embryos that Sox9 phosphorylated at S181 was localized almost exclusively in the pre-hypertrophic zone of the growth plate, an area corresponding to the major site of expression of PTH/PTHrP receptor. In contrast, no phosphorylation of Sox9 at S181 was detected in growth plates of PTH/PTHrP receptor null mutant mice. Sox9, regardless of phosphorylation state, was present in all chondrocytes of both genotypes except in hypertrophic chondrocytes. Thus, Sox9 is a target of PTHrP signaling and the PTHrP-dependent phosphorylation of SOX9 enhances its transcriptional activity. ^ In order to investigate the in vivo function of Sox9 phosphorylation by PKA, we are generating a mouse model of mutant Sox9 harboring point mutations in two PKA phosphorylation sites. Preliminary results indicated that heterozygous mice containing half amount of mutant Sox9 that can not be phosphorylated by PKA have normal skeletal phenotype and homozygous mice are being generated. ^ Lc-Maf encodes an extra ten amino acids at the carboxyl terminus of c-Maf and contains a completely different 3′ untranslated region. The interaction between SOX9 and Lc-Maf was further confirmed by co-immunoprecipitation and GST-pull down assays, which mapped the interacting domains of SOX9 to HMG DNA binding domain and that of Lc-Maf to basic leusine zipper motif. In situ hybridizations showed that RNA of Lc-Maf coexpressed with those of Sox9 and Col2a1 in areas of mesenchymal condensation during the early stages of mouse embryo development. A DNA binding site of Lc-Maf was identified at the 5′ part of a 48-bp Col2a1 enhancer element near the HMG binding site of SOX9. Lc-Maf and SOX9 synergistically activated a luciferase reporter plasmid containing a Col2al enhancer and increased the transcription of endogenous Col2a1 gene. In summary, Lc-Maf is the first identified SOX9-interating protein during chondrogenesis and may be an important activator of Col2a1 gene. ^

Analysis of Agrobacterium tumefaciens VirB6 and VirB9 of the VirB /D4 type IV secretion system

Agrobacterium tumefaciens uses the VirB/D4 type IV secretion system (T4SS) to translocate oncogenic DNA (T-DNA) and protein substrates to plant cells. Independent of VirD4, the eleven VirB proteins are also essential for elaboration of a conjugative pilus termed the T pilus. The focus of this thesis is the characterization and analysis of two VirB proteins, VirB6 and VirB9, with respect to substrate translocation and T pilus biogenesis. Observed stabilizing effects of VirB6 on other VirB subunits and results of protein-protein interaction studies suggest that VirB6 mediates assembly of the secretion machine and T pilus through interactions with VirB7 and VirB9. Topology studies support a model for VirB6 as a polytopic membrane protein with a periplasmic N terminus, a large internal periplasmic loop, five transmembrane segments, and a cytoplasmic C terminus. Topology studies and Transfer DNA immunoprecipitation (TrIP) assays identified several important VirB6 functional domains: (i) the large internal periplasmic loop mediates interaction of VirB6 with the T-DNA, (ii) the membrane spanning region carboxyl-terminal to the large periplasmic loop mediates substrate transfer from VirB6 to VirB8, and (iii) the terminal regions of VirB6 are required for substrate transfer to VirB2 and VirB9. To analyze structure-function relationships of VirB9, the phenotypic consequences of dipeptide insertion mutations were characterized. Substrate discriminating mutations were shown to selectively export the oncogenic T-DNA and VirE2 to plant cells or a mobilizable IncQ plasmid to bacterial cells. Mutations affecting VirB9 interactions with VirB7 and VirB10 were localized to the C- and N- terminal regions respectively. Additionally, “uncoupling” mutations identified in VirB11 and VirB6 that block T pilus assembly, but not substrate transfer to recipient cells, were also identified in VirB9. These results in conjunction with computer analysis establish that VirB9, like VirB6, is also composed of distinct regions or domains that contribute in various ways to secretion channel activity and T pilus assembly. Lastly, in vivo immunofluorescent studies suggest that VirB9 localizes to the outer membrane and may play a role similar to that of secretion/ushers of types II and III secretion systems to facilitate substrate translocation across this final bacterial barrier. ^

Complete Genome Sequence of Treponema paraluiscuniculi, Strain Cuniculi A: The Loss of Infectivity to Humans Is Associated with Genome Decay.

Treponema paraluiscuniculi is the causative agent of rabbit venereal spirochetosis. It is not infectious to humans, although its genome structure is very closely related to other pathogenic Treponema species including Treponema pallidum subspecies pallidum, the etiological agent of syphilis. In this study, the genome sequence of Treponema paraluiscuniculi, strain Cuniculi A, was determined by a combination of several high-throughput sequencing strategies. Whereas the overall size (1,133,390 bp), arrangement, and gene content of the Cuniculi A genome closely resembled those of the T. pallidum genome, the T. paraluiscuniculi genome contained a markedly higher number of pseudogenes and gene fragments (51). In addition to pseudogenes, 33 divergent genes were also found in the T. paraluiscuniculi genome. A set of 32 (out of 84) affected genes encoded proteins of known or predicted function in the Nichols genome. These proteins included virulence factors, gene regulators and components of DNA repair and recombination. The majority (52 or 61.9%) of the Cuniculi A pseudogenes and divergent genes were of unknown function. Our results indicate that T. paraluiscuniculi has evolved from a T. pallidum-like ancestor and adapted to a specialized host-associated niche (rabbits) during loss of infectivity to humans. The genes that are inactivated or altered in T. paraluiscuniculi are candidates for virulence factors important in the infectivity and pathogenesis of T. pallidum subspecies.

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^